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Abstract 

In this paper we investigate the dynamics of solitons occurring in the 
nonlinear Schroedinger equation when a parameter h — ► 0. 
We prove that under suitable assumptions, the the soliton approxi- 
mately follows the dynamics of a point particle, namely, the motion 
of its bary center qh(t) satisfies the equation 

q h (t)+VV(q h (t))=H h {t) 

where 

sup \H h (t)\ ->0 as h^O. 



Mathematics subject classification. 35Q55, 35Q51, 37K40, 37K45, 
47J35. 

Keywords. Soliton dynamics, Nonlinear Schroedinger Equation, or- 
bital stability, concentration phenomena, semiclassical limit. 



Contents 

1 Introduction 

1.1 Notations 



*Dipartimento di Matematica Applicata, Universita degli Studi di Pisa, Via F. Buonar- 
roti 1/c, Pisa, ITALY, e-mail: benci@dma.unipi.it, ghimenti@mail.dm.unipi.it, 
a.micheletti@dma.unipi . it 



1 



General features of NSE 8 

2.1 The first integrals of NSE Jg 

2.2 The swarm interpretation 11 

2.3 An heuristic proof 1 12 

Preliminary results |l3j 

3.1 Existence and dynamics of bary center 13| 

3.2 Concentration results TH 

The final result E 

4.1 Barycenter and concentration point '. 

4.2 Equation of the travelling soliton L' 



1 Introduction 

Roughly speaking a solitary wave is a solution of a field equation whose 
energy travels as a localized packet and which preserves this localization in 
time. 

By soliton we mean an orbitally stable solitary wave so that it has a 
particle-like behavior (for the definition of orbital stability we refer e.g. to 
PP, 0, etc.). 

In this paper we will be concerned with the dynamics of solitons relative 
to a class of nonlinear Schroedinger equations (NSE). 

Let us consider the following Cauchy problem relative to the NSE: 

i,(0,x) = ±;u(^)ei™ (2) 

where 



2 

and U : R^ —>■ R, iV > 2, is a positive, radially symmetric solution of the 
static nonlinear Schroedinger equation 

-AU + W'(U) = 2uU (3) 

with 

\\u\y = o (4) 
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Direct computations show that a solution of ([I]),© is given by 

1 



with 



„ 1 9 LO 

E = -v + i 

2 /i a -7 



(5) 



Moreover if the problem ([I]),© is well posed this is the unique solution. 

We can interprete this result saying that the barycenter q(t) of the solution 
of (EE]) defined by 

x\ip(t, x)\ 2 dx 

(6) 



satisfies the Cauchy problem 



\ip(t,x)\ dx 



q = 

9(0) = % 
q(0) = v 

The aim of this paper is to investigate what happens if the problem is 
perturbed namely to investigate the problem 



ih t = "T A * + W W '^W\ + V(x)i ' 

if) (0,x) = if h {x) 



(Ph) 



where 



<Ph{X) 



1 ( TT . x fx-qo 



en 



(7) 



and Wo is small, namely there is a constant C such that 

V{x)\w {x)\ 2 dx < Ch a ~^ 



Also we assume that 



\U + w \ 



L- 



\u\ 



L 2 



a 



We make the following assumptions: 
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(i) the problem QP/J has a unique solution 

i> G C°(R,H 2 (R N )) nC l (R,L 2 (R N )) (8) 

(sufficient conditions can be found in Kato |14| . Cazenave [Bj, Ginibre- 
Velo [HI; see Remark [2]). 

(ii) : R + — > R is a C 3 function which satisfies the following assumptions: 

W(0) = W(0) = W"(0) = (Wo) 

\W"(s)\ < c 1 \s\ q - 2 + c 2 \s\ p ~ 2 for some 2<q<p<2*. (Wi) 

4 

W(s) > -c|s| y , c > 0, 2 < z/ < 2 + — and s large (W 2 ) 
3s E R + such that W{s ) < (W 3 ) 

(iii) V : R — > R is a C 2 function which satisfies the following assumptions: 

V(x) > 0; (Vb) 



\VV{x)\ < V{x) b for |x| > R 1 >l,be (0,1); (^ 



K(x) > \x\ a for |x| > Pi > 1, a > 1. (V 2 ) 
The main result of this paper is the following theorem: 
Theorem 1. Assume that (i), (ii) and (iii) hold and that 

a > 7 (crucial assumption) 



Then the barycenter quit) of the solution of the problem (fP fc | ) satisfies the 
following Cauchy problem: 

q h (t) + VV(q h (t)) = H h (t) 

Qh(0) = qo 
q h (0) = v 

where 

sup \Hh(t)\ — > as h — > 
Let us discuss the set of our assumptions 
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Remark 2. About the assumption (i), we recall a result on the global ex- 
istence of solutions of the Cauchy problem ( T^j ) (see P, HH Ej). Assume 
dTO , ( TO and O for W. Let D(A) (resp. L>(,4 1/2 )) denote the domain 
of the selfadjoint operator A (resp. A 1 / 2 ) where 



1/21 



A = -A + V : L 2 (R V ) -> L 2 (R V ). 

If V > 0, V G C 2 and |<9 2 1/| G L°° and the initial data i/}(0,x) G £>(A 
then there exists the global solution ip of ( P/J and 

^(t, x) G C° (R, D(A 1/2 )) n C^R, F-^E^)). 
Furthermore, if ip(0,x) G -D(A) then 

ip(t, x) G C°(R, n C\R, L 2 (R W )). 

In this case, since D(A) C H 2 (R N ), (i) is satisfied. 

Remark 3. The conditions <\Wp\j and flV^j ) are assumed for simplicity; in fact 
they can be weakened as follows 

W"(0) = E 

and 

V(x) > E x . 

In fact, in the general case, the solution of the Schroedinger equation is 
modified only by a phase factor. 

Remark 4. In [2] the authors prove that if (ii) holds equation (JTJ admits 
orbitally stable solitary waves having the form (j2J). In particular the authors 
show that, under assumptions QWiD , <\W 2 \j and <\W 3 \j , for any a there exists a 
minimizer U(x) = U a (x) of the functional 



J{u) = J Q|Vm| 2 + W(u)^j dx 



on the manifold S a := {u E H 1 , ||w||ia = a}. Such a minimizer satisfies 
eq.(E]) where 2u is a Lagrange multiplier. We will call ground state solution 
a minimizer radially symmetric around the origin. We recall that by a well 
known result of Gidas, Ni, and Nirenberg [TU] , any positive solution of eq. © 
is radially symmetric around some point. 
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Remark 5. We set 

u h {x) = h^U (JL) 

where U is a ground state solution. Now we establish a relation between 
a, (3 and 7 in order to have stationary solution of ([TJ) of the form ip(t,x) = 
Uh{x)e l ~ , namely, Uh(x) is a solution of the equation 

-h 2 Au h + ^W'iWuh) = 2u h u h . 
In fact, replacing Uh by its explicit expression, we get 

- k ^(£)] + ^(£))=^-M£) 

and hence, by rescaling the variable x, 

_^2-7-2/J+«a?7(x) + W {U{x)) = 2u h h a ^U (x) . 
Thus, it is sufficient to take 

= 1 + ^ (9) 

and 

^ = ^ (10) 
to obtain the claim. In the following we always assume ([9]). 



Remark 6. The assumption ( [V2I ) is necessary if we want to identify the 
position of the soliton with the barycenter ([6]). Let us see why. Consider a 
soliton ip(x) and a perturbation 



ipd{x) = ip(x) + <p (x — d) , d G 



Even if (x) <ti ip{x), when d is very large, the "position" of ip(x) and 
the barycenter of ipd{x) are far from each other. In Lemma [23J we shall 
prove that this situation cannot occur provided that hold. In a paper in 
preparation, we give a more involved notion of barycenter of the soliton and 
we will be able to consider other situations. 

Remark 7. We will give a rough explanation of the assumption a > 7 which, 
in this approach to the problem, is crucial. In section 12.11 we will show that 
the energy Eh of a soliton iph is composed by two parts: the internal energy 
Jh and the dynamical energy G. The internal energy is a kind of binding 
energy that prevents the soliton from splitting, while the dynamical energy 
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is related to the motion and it is composed of potential and kinetic energy. 
As h — > 0, we have that (see section I2TT]) 



■h M = h Nt3 - a -~* 



and 



Then, we have that 



h a ~ 



Jh 0M 

So the assumption a — 7 > implies that, for h <C 1, G (iph) "C Jh(tph), 
namely the internal energy is bigger than the dynamical energy. This is the 
fact that guarantees the existence and the stability of the travelling soliton 
for any time. 

As far we know, this is the first paper in which there is a result of type 
Th. [T]for all times tel. However there are other results which compare the 
motion of the soliton with the solution of the equation X(t) +W(X(t)) = 
for i e [0, T] for some constant T < 00. 

Earlier results for pure power nonlinearity and bounded external poten- 
tial are in [I]. The authors have shown that if the initial data is close to 



U{^—^ L )e % h in a suitable sense then the solution iph{t, of fl^T ) satisfies for 

te[0,T\ 



^\Mt,x)\ 2 



Tjv / \A{t,x)\ 2 dx5 X (t) 
h N Jrn 



as h — > 0. 



(li; 



c 1 



Here C 1 * is the dual of C 1 and X(t) satisfies \X{t) = W(X(t)) with X(0) = 
g , A(0) = v . 

In related papers [15] and [16] there are slight generalizations of the above 
result. Using a similar approach, Marco Squassina [19] described the soliton 
dynamics in an external magnetic potential. 

In [7] and [H] the authors study the case of bounded external potential V 
respectively in L°° or confining. 

A result comparable with Theorem EH is contained in [8]. The authors 
assume the existence of a stable ground state solution with a null space non 
degeneracy condition of the equation 
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The authors define a parameter e which depends on \i and on other parame- 
ters of the problem. Under suitable assumptions they prove that there exists 
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T > such that, if the initial data ip°(x) is very close to e tp °'^ x a °^ + * 7o ?7 MO (x — 
do) the solution tp(t, x) of problem (Pi) with initial data is given by 



with < e, p = —W(a) + o(e 2 ), a = 2p + o(e 2 ) with < t < j for e 

small. 

The main differences with our result are the following. First of all we do 
not have any limitation on the time t. Also, we have an explicit estimate 
on q (which roughly speaking correspond to a). Our assumption on the 
nonlinearity W are explicit (namely flWiD , QW2D , <\W 3 \j ) and we do not require 
the null space condition which are, in general, not easy to verify. 

1.1 Notations 

In the next we will use the following notations: 

Re(z), Im(z) are the real and the imaginary part of z 
B{xq, p) = {x E : \x — xq\ < p} 



ip(t, x) = e 



,ip(t)-{x-a{t))-Hf(t) 



r)n{t){x - a(t)) +w(t) 



(13) 



B(x ,p) c =R N \B(x ,p) 
S a = {u E H 1 : ||M|| i 2 = a} 
J c h = {uEH l : J h (u)<c} 
U q (x) = U(x - q) 




\d a V(x)\ = sup 




where a E N, i 



...,i a E {1,...,N} 



I a 2 = inf J(u) = m 

u&m, Ju 2 =a 2 



2 General features of NSE 



Equation QP/J is the Euler-Lagrange equation relative to the Lagrangian 



density 



C = ihdtW - — |V^| 2 - W h (if>) - V(x) (14) 



where, in order to simplify the notation we have set 



1 
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Sometimes it is useful to write ip in polar form 

if>(t,x) = u{t,x)e iS{t ' x)/h . 



(15) 



Thus the state of the system ip is uniquely defined by the couple of variables 
(u, S). Using these variables, the action S = j Cdxdt takes the form 



S(u,S) = - 



y |Vw| 2 + W h (u) + (d t S + l - |VS| 2 + V(x) ) u 2 



dxdt 
(16) 



and equation becomes: 

- ^-Au + W' h (u) + (d t S + X - |V5| 2 + V(x)^j u = (17) 
d t (u 2 ) + V • (u 2 \7S) = (18) 
2.1 The first integrals of NSE 

Noether's theorem states that any invariance for a one-parameter group of 
the Lagrangian implies the existence of an integral of motion (see e.g. [S]). 

Now we describe the first integrals which will be relevant for this paper, 
namely the energy and the " hylenic charge" . 

Energy The energy, by definition, is the quantity which is preserved by the 
time invariance of the Lagrangian; it has the following form 



EM 

Using (TI5]) we get 



y |W| 2 + W^) + U(X)H 2 



dx 



(19) 



2 



Vu\ 2 + W h (u)J dx + J Q \VS\ 2 + V(x)J u 2 dx 



(20) 

Thus the energy has two components: the internal energy (which, 
sometimes, is also called binding energy) 



Jh(u) = 

and the dynamical energy 
G{u, S) = 



y \VU\ 2 + W h (u) ] dx 



X - |VS| 2 + V(x) ) u 2 dx 



(21) 



:22i 
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en- 



which is composed by the kinetic energy | J \ VS\ u 2 dx and the pot 
tial energy J V(x)u 2 dx. 

By our assumptions, the internal energy is bounded from below and 
the dynamical energy is positive. 

Hylenic charge Following [I] the hylenic charge, is defined as the quantity 
which is preserved by by the invariance of the Lagrangian with respect 
to the action 

if) i — ► e ie ip. 

For equation flP^D the charge is nothing else but the L 2 norm, namely: 



= J \4>\ 2 dx = J u 2 dx 

Now we study the rescaling properties of the internal energy and the L 2 
norm of a function u(x) having the form 

u(x) := h~^v (Jjp 

We have 

IHIi 2 = ^ 27 /^(|j) 2 ^ = ^- 27 /^0 2 ^ = ^- 27 IIHIi 2 . 

and 



Jh(u) 



h 2 1 
— \Vu\ 2 + T -^W{h?u)dx 



2 

h 2 ~ 2 ^ 



V x v 



x 



—W (v (±- 



dx 



^V/3+2-2 7 -2/3 



I 



using the fundamental relation (jHJ). 

Remark 8. If we choose N/3 — 2j = 0, the L 2 norm does not change by 
rescaling. This implies that the dynamical energy G, for h small, changes 
very little. 
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2.2 The swarm interpretation 

In this section we will suppose that the soliton is composed by a swarm of 
particles which follow the laws of classical dynamics given by the Hamilton- 
Jacobi equation. This interpretation will permit us to give an heuristic proof 
of the main result. 

First of all let us write NSE with the usual physical constants m and H: 

Here m has the dimension of mass and H, the Plank constant, has the di- 
mension of action. 

In this case equations (iT7|) and (ITH1) become: 

- |^Aw + l -W' h {u) + Us + ^ \VS\ 2 + V(x)j u = (23) 

Mu 2 )+V-(^)=0 (24) 

The second equation allows us to interprete the matter field to be a fluid 
composed by particles whose density is given by 

PH = U 2 

and which move in the velocity field 

VS 
v = . 

m 

So equation f[2~4"|) becomes the continuity equation: 

dtpn + V ■ (pnv) = 0. 

If 

- |^ + \w' n {u) « «, (26) 
equation (J2"3"|) can be approximated by the eikonal equation 

^5 + -^|VS| 2 + r(x) = 0. (27) 
2m 

This is the Hamilton- Jacobi equation of a particle of mass m in a potential 
field V. 



(25) 
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If we do not assume ([211 . equation (j2*7l) needs to be replaced by 

d t S+^-\\/S\ 2 + V + Q(u) = (2* 
2m 

with 



2m 

The term Q(tt) can be regarded as a field describing a sort of interaction 
between particles. 

Given a solution S(t, x) of the Hamilton- Jacobi equation, the motion of 
the particles is determined by eq.(|25l). 

2.3 An heuristic proof 

In this section we present an heuristic proof of the main result. This proof 
is not at all rigorous, but it helps to understand the underlying Physics. 

If we interpret pn = u 2 as the density of particles then 

H — pndx 



is the total number of particles. By ( l28l) . each of these particle moves as a 
classical particle of mass m and hence, we can apply to the laws of classical 
dynamics. In particular the center of mass defined in ([6]) takes the following 
form: 

= Jxmp H dx = fx PH dx 
J mp-j-cdx J pndx 

The motion of the barycenter is not affected by the interaction between 
particles (namely by the term (!28l)). but only by the external forces, namely 
by VV. Thus the global external force acting on the swarm of particles is 
given by 



VV(x)p H dx. (30) 

Thus the motion of the center of mass q follows the Newton law 

~F = Mq, (31) 

where M = f mpndx is the total mass of the swarm; thus by (|2"91) . (IHDl and 
(I3"T|) . we get 

f Wpndx f Wu 2 dx 



J pndx m J u 2 dx 
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If we assume that the u(t,x) and hence pn(t,x) is concentrated in the 
point q(t), we have that 



and so, we get 

mq(t) = -W(g(t)) . 

Notice that the equation mq(t) = — W (q(t)) is the Newtonian form of 
the Hamilton- Jacobi equation (j2j~ 



3 Preliminary results 

In this section we prove two results which are the base of the main theorem. 
They have some interest in themselves and require less assumptions than the 
final theorem. 

3.1 Existence and dynamics of barycenter 

We recall the definition of barycenter of ip 



x\tp{t,x)\ dx 

<//,(/) ^ — (32) 

\ip(t, x)\ 2 dx 



The barycenter is not well defined for all the functions ip G ff 1 (lR iV ). Thus 
we need the following result: 



Theorem 9. Let ip(t,x) be a global solution of |7Jp such that ip(t,x) G 
C(R,fl' 1 (R JV )) D C rl (R,fr -1 (R JV )) with initial data ip(0,x) such that 

\x\ \ip(0, x)\ 2 dx < +oo. 

Then the map qh{t) : R — > R , given by 13^) is well defined. 

Proof. We show that | ■ | 1,/2 |?/>(t, -)| G L 2 (M. N ) for any t, using a regular- 
ization argument. 
We set 

k £ {t) = [ e- 2£lxl \x\\ij(t,x)\ 2 dx. 
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Since ip is a solution of Q-P/J , we have 

A£(t) = / e _2e|a,| |a;|[8t^ ~ MV^] = 21m ( / ijxl^t^e -2 ^ 

= Mm ^ V^Vd^l^e- 2 ^)^ + ^Im H 2|x| VHV 2£N 



+ 



■2e a; 



= him (^J VV>VV>(Me _2eW n + him (j ^V^V(|x|e" 
= Mm (J$Vtj)- |^e- 2e|x| (l -2e|x|)J , 
so we have 

|^)| < / |^||V^|<||^,-)IIl 2 ||V^,-)IIl 2 , (33) 

then by (I3U]) we get 




x 



k £ {t) = fc £ (0) + Mm( / / ^V^-f7e- 2e|x,| (l-2£|x|)rfxc/t ) < 



'o 



< nvi^^(o,x)iii 2 + / \m,-)\\ L2 \\vip(t,-)\\ L2 dt. 



By Fatou Lemma, when e->0we get | ■ | 1/,2 |?/>(t, ■) G L 2 (R N ) for any £ > 0. 
So the map q{t) : [0, oo) — > R^ is well defined. 
□ 

Theorem 10. JTie map q h (t) : R -> R N , given by §EE) is C 1 and 



Im I hi ip(t,x)V^{t,x)dx 

V Jr n 



*M - 11,(^)111, — '- <"> 

Moreover if ip(t,x) E C(R, P 2 (R Ar ))nC 1 (lR, L 2 (R N )) then q h (t) is C 2 and 

V(x)V\ip(t,x)\ 2 dx 

*w ^ JM " p» • (35) 
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Proof. We have 



him 



\mt,x)\\i* 



We use the same regularization argument of Th. EH We set 

Klit) = [ e- 2e ^Xi\i)\ 2 dx 



and again we find in the same way that 

-him (^j ^W> ■ ^25Xi e- 2£|x| ^ 

where is the i-th vector of the canonical base of M. . So, there exists a 
constant c > such that 



< C ||V(V)IMIW(i, 



\L 2 , 



Then we have that 



A" ; '(/) - K'AO) + Mm ( j / rV r ■ < , <~- : ' - 



\x\ 



-2e\x\ 



Using the dominated convergence theorem, when e^Owe have 



Xi\ijj(t,x)\ dx = j Xi\ifj(0,x)\ dx + j him ( / ipVip ■ ei ) dt, 



so, for all i we have 

d_ 

dt 



/ Xi\ijj\ 2 dx = him I / ipVip ■ e, 

Jr n \J 



(36) 



This proves the first part of the theorem. 

( V{x)V\i){t,x)\ 2 dx 
Next we prove that q(t) = — M ; u ^ l|2 under the supplemen- 
tary assumption ip G C 1 (1R, H 1 ). 



\m, x )\\ 



L 2 
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By this assumption we have that ip(t, x)Vip(t, x) G C 1 (M, L 1 (IR Ar )). Thus 
h / lm(d t [4>(t,x)Vil>(t,x)})dx 

Jrn 



h / Im(<9^(t, x) Vip(t, x) + ip(t, x)d t Vip(t, x))d 
Jr n 











h / lm(dt^p(t, x 
Jrn 


)Vip(t,x) 


+ $(t,x)Vd t ip(t, 


x))dx 




Mt,x] 


Wh 




h / lm(d t 4>(t, x 


)W>(t,x) 


-V$(t,x)dtif>(t, 


x))dx 




||^, x) 


\\h 




2Re ( f 

\J M A 


ihd t ip(t, 


x)Vijj(t, x)dx^j 






\mt,x] 







Re 



-h 2 A^ + 2V^ + — W'(\ip\) 



Vipdx 



\mt,x)\\ 



L 2 



We have, for alH = 1, . . . , N, 

1> 



Re(/ W'M)^-d x ^ )= f d Xi WM\) = 0, 
\Jr n \V\ J Jr n 

because W(\^\) G L 1 (R JV ) and W\\^\)d x ^ G L l because ijj(t, •) G # 2 . 
In the same way we have 



Re / -A4>d Xi 4> 



Thus 



q h {t) 



a 



2Re / Vi/jVipdx) Re 



/ v(x)vm, 



x)\ 2 dx 



\m,x)\\ 2 



L 2 



\m,x)\\ 2 



(37) 



L 2 



We point out that V\tp\ G L 1 because ip is a global solution with ■?/> G -ff 2 , 

G L 2 . 
Step 3. Conclusion. 
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Let ip(t,x) G C°(R,H 2 ) n ^(R,^ 2 ). We define a function 7 A (f,x) G 
C°(R,i/ 2 ) HC^R,// 1 ) as 

7a(*,x) = / ^ A (x-0^,0^ (38) 

where </?(£) = V?(lf I) i s a positive smooth function with compact support in 
|f | < A, with / <p x (S)dt = I. 
Fixed t we have that 

-fx(t,x) -> in H 2 (R N ) for A -> 0; 

dax(t,x) -> <9^(t,x) in L 2 (IR iV ) for A -> 0, 

and the convergence is uniform for every compact set in R. 

Furthermore, using that ^ is a global solution in C°(R, H 2 ) n C 1 (1R, L 2 ) 
we have that V(x)^(*, x) € C°(R,L 2 ) and that, fixed t 

V(x) lx (t,x) -> \/(x)^(t,a;) in L 2 (M Ar ) for A -> 0; 

again, the convergence is uniform for every compact set in R. 
We have that 7 A (t,x) solve the following differential equation 

= + ^'(I7a(M)|)^|^ + ^(x) 7 a(«,x) + r( 7A (M)) 

where, r(7 A (t,x)) — > in L 2 , for all t, as A — > 0, uniformly on every compact 
set in R. Thus we have, proceeding as in Step 3, 



eft 



Im(7 A (t,x)7 A (t,a;)) = 

V(x) V|7 A (t, x)| 2 + 2Re(r( 7A (t, x)) V 7 (t, x))da:. (39) 



/ 



/R^ 

We have 



/ Im(7 A (t,x) 7A (t,a:)) = / Im(7 A (0, x) 7a (0, x))+ 

+ 11 V(x)\7\-f X (s,x)\ 2 + 2Re(ri-f X (s,x))\7^(s,x))dxds (40) 
Jo Jr n 
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and, for all s, 

V(x)V\^x(s,x)\ 2 + 2Re(r( lx (s,x))V^(s,x))dx -> / V(x)V\ip(s, x)\ 2 

Jr n 

as A — > 0. Finally, 

y V(x)V|7 A (s,x)| 2 + 2Re(r(7 A (s,x))V7(s,x))c/x < 

<||^)7A( S r)IUHl7A( S r)ll^ + ll(r(7A(«r))IUHl7A( S ,-)ll^- (41) 

And because V(x) r yx(s, •) — > V(a;)^(s, •) in L 2 , r(7 A (s, •) -> in L 2 and 
7a(s, •) — > ip(s, •) in if 1 uniformly in s on every compact set, we have that 
for some constant C 

sup / y(x)V|7 A (s,x)| 2 + 2Re(r(7 A (s,x))V7(s,x))rfx < C. (42) 
se[o,t] Jr^ 

By (01]) we get 

/ V r (x)V|7 A (s, x)| 2 c?X(is + / / 2Re(r(7 A (s, x))V7(s, x))dxds — > 
o it* Jo Jr^ 

/ I V{x)V\ip{s,x)\ 2 dxds. (43) 
Jo Jr^ 

Furthermore we know 



Im( 7A (t,x)7 A (t,x)) -> / Im(^(t,x)^(t,x)); (44) 

Im(7 A (0,x)7 A (0,a;)) -> / lm(^(0, x)^(0, x)). (45) 

Jm N 

At this point by flUD, (PJ, (JSj) and (US]) we get 
In^^t, x)if)(t, x)) = 

= [ lm(^(0,a;)^(0,a;))+ / / V(a;)V|^(s, x)| 2 cfocfe (46) 

JR^ JO JR^ 

that concludes the proof. 

We have the following corollary 
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Corollary 11. Assume (ViP and the assumptions of the previous theorem; 
then 

[ W(x)\ip(t,x)\ 2 dx 
= ~ n:u m,2 • (47) 



\m, X )\\ 



L' 2 



Proof. By §Vj]) , we have that 



< / V\x)\ij(t,x)\ 2 dx 



VV{x)\ip{t,x)\ 2 dx 



< Ci / V(x)\ip(t,x)\ 2 dx 

< C 2 G < +00 

where G is the dynamical energy (1221) . Thus, we can integrate by parts and 
we have that 



V(x)V\ij(t,x)\ 2 dx 



VV(x)\ij(t,x)\ 2 dx 



□ 



Remark 12. If we use the polar form (fT5|) . 
meaningful form respectively: 



and (I47p take the more 



Qh(t) 



Qh(t) 



VS u dx 



u dx 



W{x)u 2 dx 



u dx 



They can be interpreted as follows: qh{t) is the average momentum (remem- 
ber ( 1251) and that m = 1); %(t) equals the average force, since F = —W 
(see dSUD). 



3.2 Concentration results 



In this section we prove a concentration property of the solution of flP^D with 
initial data (JTj); more exactly, we prove that for any time i el, this solution 
is a "bump" of radius less than some constant " R" . In order to prove this 
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result, it is sufficient to assume that problem QP/J admits global solutions 
ip(t,x) G C(R, if 1 (IR Ar )) which satisfy the conservation of the energy and of 
the L 2 norm, namely it is not necessary to assume the regularity (E])- 

For some nonlinearities W, it is possible that the ground state solution is 
not unique. In any case we have the following result: 

Proposition 13. Let U be a ground state solution of Then, for \x\ > 1 
and N > 2 

C 

U{x) < N _, 

\X\ 2 

where C is a constant which does not depend on U . 

Proof. By a well known inequality due to Strauss [2U] we have 



< U(x) < c j^jj^ 1 for \x\ > a N a.e. (4f 



\X\ 2 



where Cn and depend only on the dimension N. Moreover there exists a 
constant C m , such that < C m for any U minimizer of inf J(u) = m. 

ueScr 

In fact we have the following inequality 



\\u\\ L » < K\\u\\ L2 2 v \\Vu\\l 2 " (49) 
for some constant b v . Then, by (jlj) 

||t/||^<6>"(^)||VC/||"r A '. (50) 



By assumption JW2D and by ( J50l) we have 

1 



m = J{U)> -\VU\ 2 - Cl U 2 - c 2 U u dx> 



> ^ / |Vt/| 2 -c 3 / \VU\ 2 ) - Cl a 2 



1 

2 „ 

for some constant C3. If < — N < 2, namely 2 < v < 2 + 4, we have 
the claim. 
□ 

Lemma 14. For any e > 0, there exists an R = R(e) and a 5 = 5(e) such 
that, for any u G J m+5 n S a , we can find a point q = q{u) G M. N such that 



a 2 



u 2 (x)dx < e. (51) 

r \B(q,R) 



20 



Proof. Firstly we prove that for any e > 0, there exists a 5 such that, for 
any u G J m+S f] S a , we can find a point q = q{u) G K.^ and a radial ground 
state solution U such that 

|Kx)-C/(x-g)|| H i <£■ (52) 

We argue by contradiction. Suppose that there exists an e > and a 
sequence {w n }n such that ||w n ||z, 2 = cr, J(it n ) — > m and, for any q G and 
for each U ground state solution it holds 

s< \\u(x) — U(x — (53) 

By the Ekeland principle we can assume that {u n } is a Palais Smale 
sequence for J on S a , that is, there exists {A n } such that 

— Au n + W'(u n ) — \ n u n — ¥ as n — > oo. (54) 

By |21 Proposition 11] up to a subsequence we have that A n — > A < 0. So 
we get 

-Aw„ + W\u n ) - Aw n -> 0; (55) 
J(u n ) — A / — > m — Act 2 . (56) 



As a consequence of the Concentration Compactness principle [17, 18J, we 
can describe the behavior of this P.S. sequence. We use the Splitting Lemma 
(see [21], [2]) and we get 



k 

U n = 

k 



U j {x-ql)+w n with w n -> in H 1 (57) 

j=i 

k 

^ = Ell^-gjl)lli» (58) 

it 

^J(^(x-g^))=m = / (T2 (59) 



where are solutions of —Am + W'{u) = Xu and q n G 

Here 
we have 



DiV 



Here J p 2 = min|| u ||2 2=p2 J(u). We recall (see [T7]) that for any G (0, p) 



I p 2 < I U 2 + J p 2„ M 2. (60) 
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We verify that in (j571) - (j591) it is k — 1. We assume k = 2. Suppose that 
HC/ 1 1 1 2,2 = /i 2 < a 2 . Then, by ([50]) . we have a contradiction because 

1,2 < V + 1,2^2 < J{U l ) + J(U 2 ) = 1,2. (61) 

For the case k > 2 we argue analogously. 
Thus we have, up to subsequence, 

u n (x) = U(x - q n ) +w n w n -> in (62) 

for some [/ radial ground state solution, and ( |62|) contradicts ( |53l) . 

At this point, given e, there exist a point q = q(u) G M. N and a radial 
ground state solution U such that 

u(x) — U(x — q) + w and | |iu| |#i < Ce. (63) 

Now, we choose R such that 

y £/ 2 (x)cfe < (64) 

R JV \B(0,R) 

for all C/ radial ground state solutions. This is possible because, if U is a 
radial minimizer of J(u) on S a , then, as showed in following Remark | 



U(x) < C{m N l N ' ) for |x| » 1, N > 2, (65) 



X 2 

the constant C (m, N) depending only on the dimension N and on m = 
inf J(u). 

We get 

— / u 2 {x)dx <C — r / £/ 2 (a; — g)<ix H — - / w 2 + 2wUdx 

<y J v J o J 

B(q,R)C B(q,R)C B(q,R)C 

= J U 2 (x)dx + -^ J w 2 + 2wUdx. (66) 

B(Q,R) C B(q,R) c 

By (l6"3"j) . (!B^|) . (I6"B"]) we get the claim. We notice also the R does not depend 
on u, U. 
□ 



We can describe now the concentration properties of the solution of (\Ph\j - 
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Lemma 15. For any e > 0, there exists a 5 = 5(e) and a R = R(e) such 
that for any ip(t, x) solution of /jP h \) with \h^ip{t, hPx)\ G J m+S n S a for all 
t there exists a qu{t) G ^ N for which 



1 



a 2 h Nf3-2j 



\ip(t, x)\ 2 dx < e. 



(67) 



Here qh(t) depends on e and ip(t,x). 

Proof. Fixed h and t, we set v(£) = \h~ lr il)(t, hr£)\. So we have 
m < J(v) < m + 5 and | M|l 2 — °~- 

So, by Lemma [HI we have that there exist an R > and a q = q(v) such 
that 

1 



6 > — 
a 2 



(68) 



WB(q,R) 



By a change of variable we obtain 



£ > — 
a 2 



i 



a 



2^/3-27 



(t,x)\ dx, 



(69) 



where quit) depends on e, h, t and ip, while R depends only by e. 
□ 

We give now some results about the concentration property of the solu- 
tions ip(t, x) of the problem (\Ph\j - Given K > q G R , h > we call 

ip(0,x) = u h (0,x)ei Sh< -°' x) 
with u h (0,x) = h-~< [(U + w) (2^)] 



U is a radial ground state solution 

| \U + w\\ L 2 = | \U\ \ L 2 = a and | \w\ \ H i < Kh a ^ 

||V5ft(0,x)|| L <- < K for all h 

j RN V{x)u 2 h {0,x)dx < Kh N ?- 2a . 
the set of admissible initial data. 



(70) 
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Remark 16. The condition ||w||^i < Kh 01 " 1 can be weakened. Indeed in 
the proof of the theorem we need J(U+w) < m + Kh a ~' y , which is implied by 
I I'M In 1 < Kh a ~' 1 . We prefer to refer to the strongest but simpler hypotheses 
to simplify the statement of the main theorem. 

Remark 17. In TheoremfHwe assume 5^(0, x) = v-x which is more stronger 
than ||V5/i(0, sc) 1 1 2,00 < K to simplify the statement and for a better physical 
interpretation. 

Finally, we can prove the main result of this section. 

Theorem 18. Assume V G Lg, and fV^]). Fix K > 0, q eR N . Let a > 7. 

For all e > 0, there exists R > and h > such that, for any ip(t,x) 
solution of with initial data ip(0, x) G B h ' q with h < h , and for any t, 
there exists qh(t) G M. N for which 



\m,x)\\ 



L-' 



\ip(t, x)\ 2 dx < e. 



(71) 



R N ^B(q h (t),RhP) 



Here quit) depends on tp(t,x) 



Proof. By the conservation law, the energy Eh(ip(t, x)) is constant with 
respect to t. Then we have 



E h {if>(t,x)) 



E h (ip(0,x)) 
J h (u h (0,x)) + 



u 2 h (0,x) 



\VS h (0,x) 



+ V(x) 



dx 



K 



< J h (u h (0, x)) + jo- 2 h N ^ + Kh N ^ 

= h Nf3 - a ~~<J(U + w) + h N P-^C 

where C is a suitable constant. Now, by rescaling, and using that i/j(0,x) G 
Bh' q , we obtain 

E h (iP(t,x)) = h N ^~ a J(U + w) + Ch N ^ 

< h N ^- a {m + Kh a -^) + Ch Np -^ (72) 
= h N ^~ a (m + Kh a ~~< + Ch a ^) = h N/3 -^~ a (m + h a -^d) 

where C\ is a suitable constant. Thus 



Jh(u hn (t,x)) = E h (ip(t,x)) - G h (ip(t,x) 
E h {il)[t,x)) 



+ V(x) 



Uh(t, x) 2 dx 



< h N0 - y - a (m + h a -' Y C 1 ) 



(73) 
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because V > 0. By rescaling the inequality (1731) we get 

J(h^u h {t, h p x)) <m + h a ^C x (74) 

So, if a > 7, for /i small we can apply Lemma [TBI and we get the claim. 
□ 



4 The final result 

4.1 Bary center and concentration point 

In this paragraph, we estimate the distance between the concentration point 
and the barycenter of a solution ip{t,x) for a potential satisfying hypothesis 



m and (EJ). 

Hereafter, fixed K > 0, we assume that ip(t,x) is a global solution of 
the Schrodinger equation M, ip(t, x) <G C(R, if 1 ) nC^l, if -1 ), with initial 
data ip{0,x) G .Bf' 9 with given by ([TOD - 

Lemma 19. There exists a constant L > swc/i t/iat 

Proof. At first we notice that H/r^w^t, ^)||| 2 = H/r^u^O, h p x)\\j^ = 
\\U + w\\ 2 L2 = a 2 . Thus 

J h (u h (t, x)) = h N/3 -^J(h^u h (t, h p x)) > h Np ~^m. (75) 

By (1721) . there exist a constant L such that 

E h (ip(t, x)) < h N ^~ a m + Lh N/3 " 2 \ (76) 

Finally, 

f |VS1 2 

\/(a;)M^(t,x)(ix = E h (ip(t,x)) - J h {u h {t,x)) -/ — —u 2 h (t,x)dx 

Jrn l 

< E h {ip{t, x)) - J h (u h (t, x)) 

< h N ^~ a m + Lh N ^ - h Nfi -^m = Lh N ^ 

that concludes the proof. 
□ 
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Remark 20. By Lemma [T9l we get, for any R 2 > R± (R\ given in (V2)) and 
for any i Gl the following inequality 



L > 



> 



1 



|x|>i?/2 



j^ct— 1 



(77) 



|x|>-R 2 



1 37 1 ^ i ^ x^jdx 



\x\>R 2 



Lemma 21. There exists a constant Ki such that 

\q h {t)\ <#i /or* 6 1. 
By Lemma [T9l and Remark |2"U1 we have that 



< I \x\u 2 h (t,x) + / \x\ul(t,x) 

'|a;|>Ri V|x|<_Ri 



< R, [ ulit,x)dx + ^ z - 1 h N ^. 
So, using the definition of qhif) we have 



lfc(*)l<*i+j^5 = *i 



(7S 



for some f^i > 0. 



Remark 22. By the inequality 0771) in Remark [201 we have also that, for 
any R 2 > R\, 



u^it, x)dx 



\x\>R 2 



< 



u1(t, x)dx 



a 2 Rl 



for all t e 



Hereafter, we always choose R 2 large enough to have 



L 1 

< 



a 2 R a 2 ' 2" 



(79) 



Now we show the boundedness of the concentration point qh{t) defined 
in Lemma (fTo"]) . 

Lemma 23. Given < e < 1/2, and R 2 as in the previous remark. 
We get 



26 



1. sup |<?/i(t)| < R 2 + R(£)hP < R 2 + 1, for all h < h and 5 < 5 small 
enough. 

2. sup teR \ q h (t) - qh(t)\ < ^k=i + 3-R 3 £ + R{e)hP , for any R 3 > R 2 , and 



for all h small enough. 



Proof. 

Stepl. We prove the boundedness of the concentration point qhif). 

By the Theorem [TBI with e < 1/2, and by Remark [2*2*1 it is obvious that 
the ball B(qh(t), R{e)h 13 ) is not contained in the set ~R N s B(0, R2), and we 
have 

B(q h (t), R(e)h^) c 5(0, R 2 + 2R{e)h p ). (80) 

Because R(s) does not depend on h, we can assume h so small that 2R{e)h li < 
1. Then 



|&(f)| < R 2 + 2R(e)h? <R 2 + 1; 
B(q h {t),R{e)h l3 )GB(0,R 2 + l); 



(81) 
(82) 



This concludes the proof of the first claim. 

Step 2. We estimate the difference between the barycenter and the concen- 
tration point. 
We have 



kh(t) -qh(t)\ 



(x - q h (t))u 2 h (t,x)dx 



and we split the integral in three parts, with R 3 > R 2 : 



h 



f (x - q h (t))u 2 h (t,x)dx 

R N ^B(0,R 3 ) 



J uf l (t,x)dx 



f (x - q h (t))u 2 h (t,x)dx 

A 2 



J u\(t,x)dx 



f (x - q h (t))u 2 h (t,x)dx 
a 3 



where A 2 = B(0,R 3 ) \ B(q h (t),R{e)h p ); 



J uf l (t,x)dx 



where A 3 = B{0,R 3 ) n B(q h (t), R{e)h p ). 
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It's trivial that -Z3 < -R(e)/?/ 3 . By Lemma I2H1 and by Theorem [T81 we have 

h < [2Rs + l]e- (84) 
By Step 1 and Remark [221 we have 

Wg(«,»)>fa < + L < _2L 
J RN uf l {t,x)dx a 2 R% a 2 Rl 



Also, by Remark 



f RN u 2 h (t,x)dx a 2 R 



< -i^T- ( 86 ) 



3 



hence 

^ ^ T^P (87) 



Concluding, we have that 



M*) - I < -TE^T + 3R z £ + R ( £ ) h ^ ( 88 ) 

for all t G K. 

□ 

We notice that Ri, R 2 and i? 3 defined in this section do not depend on e. 



4.2 Equation of the travelling soliton 

We prove that the barycenter dynamics is approximatively that of a point 
particle moving under the effect of an external potential V(x). 



Theorem 24. Assume that V satisfies (EP;(EP? (EP- Given K > 0, q E 
R N , let ip(t, x) E C(R, H 2 ) n C^M, i? 1 ) 6e a #/o&a/ solution of equation jffl) , 
wii/i initial data in B^ ,q , h < h . Then we have 

q h (t) + W(g fc (*)) = (89) 

wii/i 11^(^)11^00 goes to zero when h goes to zero. 

Proof. We know by Theorem [9j that 

Ln VV(a;)u?(i, x)dx , . 

4h(t) + jRN r \' h \/ = (90) 

Jm. N u h\ l i x ) ax 
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Hence we have to estimate the function 

S^[VV{x)-VV{q h {t))\u 2 h {t,x)dx 



H h (t) = [VV(q h (t))-VV(q h (t))\ 
We set 



M = max \d a V(r) 

a = 1, 2 
\t\ < A'i + R 2 + 1 



(91) 
(92) 



where K\ is defined in Lemma EH and i?2 is defined in Remark [22j 
By Lemma |2~T1 and Lemma I2"3l we get 



|w(&(*))- w(<fo(*))| 



< max 

i, j = 1, . . . , N 

\t\ < ifi + R 2 + 1 



d 2 V(r) 



\q h {t) - q h {t)\ < 



< M 



3L 



o 2 R 



+ 3R 3 e + R{e)h p 



(93) 



for any R 3 > R 2 . 
To estimate 



[VV{x) -VV{q h {t))]u 2 h {t,x)dx 



u\(t, x)dx 



we split the integral three parts. 



/ 



(t),R{e)h(>) 



\VV{x) - VV{q h {t))\u 2 h {t,x)dx 



M JV \ J B(g h (t),fi(£)/i' 3 ) 



u1(t, x)dx 
\W(x)\u 2 h (t,x)dx 



u1(t, x)dx 



^ N ^B{q h {t),R{e)hP) 



\VV{q h {t))\u 2 h {t,x)dx 



u 2 (t, x)dx 



By the Theorem [T8l and by Lemma [231 we have L3 < Me. 
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By (|82|) and Lemma [231 we have 



B(q h (t),BHe)hP) 



max 

i,j=l,...,N 

\t\ < A'i + R 2 + 1 



d 2 V{r) 



R{e)h^u\{t, x)dx 



< MR(e)h p . 



u1(t, x)dx 



Now, using hypothesis ( JVTT ) ? equation f[g2l . Theorem [TBI and Remark |2"U[ we 
have 



I I W( -"IMV)II" 



< 



|Vy (a; 



( u 2 h {t,x) \ b f u 2 h (t,x) 



VIM*>-)III 2 ; VIMv)il| 2 



1-6 



< 



< 



< 



B(0,R 2 +1) C 



\VV(x)\ 



lk(V)ll 2 



L 2 



6 r 



B(0,_R 2 +1) C 



IK(V)II 2 



1-6 



|w( 



.r 



uj[(t,a;) 



IMV)II 2 



-dx 



B(0,R 2 +1) C 



L 2 



.1-6 



n 6 



IMV)II 2 



L 2 



2L 



a- 1 



.1-6 



(94) 



where b G (0,1) is defined in flT^] ). Furthermore, again by Theorem [TBI we 
have 

/ I W ( x )I m I^TT^T! dx<Me. (95) 

B(0,i? 2 +l)\-B(q h (t),R(e)/i) 



IMv)ii 2 



L 2 



So, by fSD and ([95]), 

L 2 < Me + 

Concluding 

\J RN [WV(x) -VV(q h (t))]<(t^)dx\ 



.1-6 



< 2Me + 



L 



(96) 



s l - b + MR(e)h p . (97) 
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Finally, by (J23j) and ($7$ we have 

b 

£ 1 - b + M(2 + 3R 3 )e + 2MR{e)h> 3 . (98) 

At this point we can have sup \Hh{t)\ arbitrarily small choosing firstly -R3 

t 

sufficiently large, secondly e sufficiently small, and finally h small enough. 
□ 

Proof of Theorem [H By Theorem [2H we get immediately the proof of 
Theorem [H 
□ 
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